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Abstract 

We have formulated new quantum model of the QCD vacuum using the 
effective potential approach for composite operators. It is based on the ex- 
istence and importance of such kind of the nonperturbative, topologically 
nontrivial excitations of gluon field configurations, which can be effectively 
correctly described by the (?~^-type behaviour of the full gluon propagator in 
the deep infrared domain. The ultraviolet part of the full gluon propagator 
was approximated by the asymptotic freedom to-leading order perturbative 
logarithm term of the running coupling constant. Despite the vacuum en- 
ergy density remains badly divergent, we have formulated a method how to 
establish a finite (in the ultraviolet limit) relation between the two scale pa- 
rameters of our model. We have expressed the asymptotic scale parameter as 
pure number times the nonperturbative scale, which is inevitably contained 
in any realistic Ansatz for the full gluon propagator. 

PACS numbers: 11.15.Tk, 12.38.Lg 
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I. INTRODUCTION 



The modern theory of strong interactions, quantum chromodynamics (QCD) [1] is in- 
trinsically ultraviolet (UV) divergent theory because of asymptotic freedom [2,3] which is 
due to self-interaction of massless gluons only (non-Abehan character of QCD). It implies 
that any exact solution to QCD should become asymptotically free in the deep UV limit. 
This means that the perturbative UV divergences will be inevitably encountered in order to 
calculate any physical observable in this theory. Thus the existence of the UV divergences 
in QCD does not depend on the method of solution because behind these divergences is a 
physical phenomenon - asymptotic freedom. In other words, there is no way to escape from 
the UV divergences in QCD by implementing some sophisticated (other than perturbation 
theory) method of solution. Since divergences are physical (not artificial due to a method of 
solution), a mass scale parameter appears as a result of the UV renormalization program in 
order to render the theory finite in the UV limit. In QCD quark masses are not observable, 
so the UV divergences are to be absorbed into the scale parameter which precisely becomes 
finite in the UV limit. Instead of the dimensionless fixed coupling constant (like in QED, 
see a brief discussion below) we are left with the above-mentioned dimensional asymptotic 
scale parameter, IVqcd and therefore dimensionless coupling eventually becomes running. 
Numerically it is a few hundred MeV only, so it can not survive in the UV (perturbative) 
limit, i.e., it can not be determined by the theory in the UV region. It should come from 
the infrared (IR) nonpcrturbative region. This clearly shows that this scale is completely 
nonperturbative by origin and this conclusion is strongly supported by the phenomenon of 
"dimensional transmutation" [4]. 

Let us emphasize that the existence of the asymptotic scale parameter, Aqcd, rneans 
that the full gluon propagator depends in general on this scale parameter, i.e., its behavior 
in the deep UV limit deviates from the free perturbative one (the above-mentioned running 
ccfcctive coupling constant). In other words, the perturbative dynamics even in massless 
QCD is nontrivial (scale violation, asymptotic freedom). Moreover, let us note in advance, 
that due to the origin of the QCD scale parameter from the IR region, the asymptotics of the 
full gluon propagator in the deep IR limit should deviate from the free one as well. In other 
words, in the IR limit the QCD coupling is also running with its, generally speaking, own 
scale parameter (see below). In this connection, let us note that if QCD itself is confining 
theory, then some characteristic scale (different, in general, from asymptotic scale parameter, 
Aqcd) is very likely to exist. 

All this is in strict contradiction with quantum electrodynamics (QED), the theory of 
electromagnetic interactions. In QED the UV divergences are result of the method of so- 
lution, the perturbation theory. Any exact solution to QED, in principle, should not be 
asymptotically free. That is why neither the perturbative UV divergences nor the dimen- 
sional scale parameter necessarily should appear in QED. So the UV divergences in QED are 
not physical, i.e., there is no some nontrivial physics behind these divergences (unlike QCD, 
where asymptotic freedom is behind the UV divergences). However, mathematically the UV 
renormalization program (in order to render the theory finite) is, of course, the same as in 
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QCD[| but unlike QCD, all the UV divergences in QED can be absorbed into the electron 
charge and mass terms. So there is no running coupling constant and the photon propagator 
remains free one in the deep UV limit (no deviation from the perturbative free behavior in 
this limit). There is no doubt, that the drastic difference between these two theories is due 
to Abelian character of QED (no direct interaction between massless photons). 

If the structure of QCD in the UV perturbative region is well understood (because of 
asymptotic freedom, as it was briefly discussed above), little is known about its structure 
in the nonperturbative IR region. The main reason of this is, of course, the phenomenon 
of confinement which remains the most important unsolved yet problem in QCD. It is 
still remains to understand whether confinement (in particular quark confinement) is some 
additional external requirement in order to complete QCD to be a self-consistent theory of 
strong interactions or QCD itself is able to explain confinement and other nonperturbative 
phenomena without involving some extra degrees of freedom. At the very beginning of 
QCD, it was expressed an idea [5,6], that the quantum excitations of the IR degrees of 
freedom because of self-interaction of massless gluons only made it possible to understand 
confinement as well as other nonperturbative effects, for example such as dynamical chiral 
symmetry breakdown (DCSB). In other words, the importance of the IR structure of the 
true QCD vacuum was emphasized as well as its relevance to quark confinement, DCSB and 
other way around. 

In accordance with these ideas, let us suppose that QCD is intrinsically IR divergent 
theory as well. Again the physical reason of these divergences is self-interaction of massless 
gluons only. Asymptotic freedom is behind the UV perturbative divergences, while the non- 
perturbative IR divergences themselves are behind confinement and other nonperturbative 
effects in QCD. Thus in our picture they are more fundamental than the perturbative ones. 
How to take them into account? Of course, by the corresponding behavior of the full gluon 
propagator in the deep IR domain since there is no hope for an exact solution(s). In this con- 
nection, let us remind the reader that the full dynamical information of any quantum gauge 
field theory such as QCD is contained in the corresponding quantum equations of motion, 
the so-called Schwinger-Dyson (SD) equations for lower (propagators) and higher (vertices 
and kernels) Green's functions. These equations should be also complemented by the corre- 
sponding Slavnov- Taylor (ST) identities which in general relate the above mentioned lower 
and higher Green's functions to each other [1,7,8]. These identities are consequences of the 
exact gauge invariance and therefore "are exact constraints on any solution to QCD" [1]. 
Precisely this system of equations can serve as an adequate and effective tool for the non- 
perturbative approach to QCD [1,9,10]. Among the above-mentioned Green's functions, the 
two-point Green's function describing the full gluon propagator 



has a central place [1,7-11]. Here ^ is a gauge fixing parameter (^ = 0, Landau gauge) 
and T^y{q) = g^y — q^qi,/q^ = g^u — L^yi^q). Evidently, its free perturbative (tree level) 
counterpart is obtained by simply setting the full gluon form factor c?(— g^,^) = 1 in Eq. 



Though historically the renormalization of QCD was done following the renormaUzation of QED. 




(1.1) 
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(1.1). In particular, the solutions to the above-mentioned SD equation for the full gluon 
propagator (1.1), are supposed to reflect the quantum structure of the QCD ground state. It 
is a highly nonlinear integral equation containing many different propagators, vertices and 
kernels [1,7-11]. For this reason it may have many (several) different exact solutions with 
different asymptotics in the deep IR limit (the UV limit because of asymptotic freedom is 
uniquely determined), describing thus many different types of quantum excitations of gluon 
field configurations in the QCD vacuum. Evidently, not all of them can refiect the real 
structure of the QCD vacuum. Let us emphsize that we even do not know the complete set 
of boundary conditions in order to attempt to uniquely fix theory. 

The deep IR asymptotics of the full gluon propagator can be generally classified into the 
three different types: 1) the IR enhanced (IRE) or IR singular (IRS) (for references see the 
next section below), 2) the IR finite (IRF) [12] and 3) the IR vanishing (IRV) ones [13]. As 
it was emphasized above, any deviation in the behavior of the full gluon propagator in the 
IR domain from the free one automatically assumes its dependence on a scale parameter 
(at least one) responsible for the nonperturbative dynamics, say, A^vp- This is very similar 
to asymtotic freedom which requires asymptotic scale parameter, Aqcd associated with 
nontrivial perturbative dynamics (scale violation). Thus, any realistic Ansatz (since exact 
solution(s) is(are) not known) for the full gluon propagator in QCD will contain at least the 
two above-mentioned scale parameters. In our approach, the nonperturbative scale is more 
fundamental than the asymptotic one, and if the above-mentioned Ansatz is realistic one, 
than, in principle, it should give rise to the relation between them, namely the asymptotic 
scale as pure number times the nonperturbative scale. The main purpose of this article is 
precisely to establish such relation within our approach to nonperturbative QCD and its 
ground state. 

However, a few general remarks are in order. The nonperturbative QCD vacuum is a 
very complicated medium and its dynamical and topological complexity [1,14-17] means that 
its structure can be organized at various levels (classical and quantum) and it can contain 
many different components and ingredients which contribute to the truly nonperturbative 
vacuum energy density , one of the main characteristics of the QCD ground state. The clas- 
sical component of the vacuum energy density is determined by the topologically nontrivial, 
instanton-type ffuctuations of gluon fields configurations in the QCD vacuum which are non- 
perturbative weak-coupling limit solutions to the classical equation of motion in Euclidean 
space. Numerically this contribution is estimated by the random and interacting instanton 
liquid models (RILM and IILM) of the QCD vacuum (see Ref. [18] and references therein). 

Here we are going to discuss the quantum part of the vacuum energy density which is 
determined by the effective potential approach for composite operators [19,20]. It allows us 
to investigate the nonperturbative QCD vacuum, since in the absence of external sources 
the effective potential is nothing but the vacuum energy density. It gives the vacuum energy 
density in the form of the loop expansion where the number of the vacuum loops (consisting 
of the confining quarks and nonperturbative gluons properly regularized with the help of 
ghosts) is equal to the power of the Plank constant, k. 

It is well known, however, that the vacuum energy density is generally badly divergent in 
quantum field theory, in particular QCD (see, for example, the discussion given by Shifman 
in Ref. [14]). The main problem thus was how to extract the truly nonperturbative vacuum 
energy density which should be finite, automatically negative without imaginary part (stable 
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vacuum). This problem was addressed and resolved in our recent publications in Ref. [21] 
for covariant gauge QCD and in Ref. [22] for axial gauge QCD. Though the truly nonper- 
turbative vacuum energy density is important in its own right (it is the bag constant, apart 
from the sign, by definition and precisely it is related to the gluon and quark condensates, 
topological susceptibihty, etc. (see Refs. [23-26] and references therein)), here we are going 
to investigate the general structure of the Yang-Mills (YM) vacuum within our approach. 

In the next section II we will formulate our main dynamical assumption which lies in 
the heart of our approach to nonperturbative QCD in general and to its ground state in 
particular. Then in section III we will apply to the above-mentioned approach for composite 
operators in order to analyze the vacuum structure. In section IV we present our main result 
while section V is devoted to discussion. Our conclusions are given in section VI. 



II. THE IRE GLUON PROPAGATOR. ZME QUANTUM MODEL 

Today there arc no doubts left that the dynamical mechanisms of the important non- 
perturbative quantum phenomena such as quark confinement and dynamical (or equivalently 
spontaneous) chiral symmetry breaking (DCSB) are closely related to the complicated topo- 
logically nontrivial structure of the QCD vacuum [1,14-17]. On the other hand, it also 
becomes clear that the nonperturbative IR dynamical singularities, closely related to the 
nontrivial vacuum structure, play an important role in the large distance behaviour of QCD 
[5,6]. For this reason, any correct nonperturbative model of quark confinement and DCSB 
necessarily turns out to be a model of the true QCD vacuum and the other way around. 

Our model of the true QCD ground state is based on the existence and importance of 
such kind of the nonperturbative, quantum excitations of the gluon field configurations (due 
to self-interaction of massless gluons only, i.e., without explicit involving some extra degrees 
of freedom) which can be effectively correctly described by the behavior of the full gluon 
propagator in the deep IR domain (at small q^) [27,28] . Thus our main Ansatz for the full 
gluon form factor in Eq. (1.1) in the whole momentum range [0, oo) is 

Evidently, the UV piece is nothing but the running coupling constant g'^{—q'^), which is 
the leading order solution to the renormalization group equations [1] and c = (167r^/ll) 
is the inverse of the first coefficient of the /3-function in the case of YM gluon fields (pure 
gluodynamics) . In the above-mentioned Ref. [1] it is denoted as 2 /bo. The flavorless QCD 
asymptotic scale parameter is denoted as Aym- In the weak coupling limit {—q^ — > oo) 
precisely the second term becomes dominant (asymptotic freedom). 

On the other hand, we obtain the generally accepted form of the IR singular asymptotics 
for the full gluon propagator (for some references see below) 

D.Aq) - ilT', - 0, (2.2) 

which may be equivalently referred to as the strong coupling regime [1]. It describes zero 
momentum modes enhancement (ZMME) dynamical effect in QCD at large distances. In our 
previous publications and in future we prefer to use simply ZME (zero modes enhancement) 
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since we work in momentum space always. This is our primary dynamical assTimption in 
this section. The main problem with this strong singularity is its correct treatment by the 
dimensional regularization method [29] within the distribution theory [30], which precisely 
was one of highhghts of our previous pubheations [27,28] (see also Refs. [31,32]). There 
exist many arguments in favor of this behaviour. Let us remind the reader a few of them 
which arc the most important in our opinion. 

a) . Precisely this strongly singular behaviour of the full gluon propagator in the IR 
domain leads to the area law for static quarks (indicative of confinement) within the Wilson 
loop approach [33]. 

b) . The cluster property of the Wightman functions in QCD fails and this allows such 
singular behavior like (2.2) for the full gluon propagator in the deep IR domain [34]. 

c) . After the pioneering papers of Mandelstam in the covariant (Landau) gauge [35] 
and Baker, Ball and Zachariasen in the axial gauge [36], the consistency of the singular 
asymptotics (2.2) with direct solution to the SD equation for the full gluon propagator in 
the IR domain was repeatedly confirmed (see for example Refs. [9,10, 37,38] and references 
therein) . 

d) . Moreover, let us underline that without this component in the decomposition of the 
full gluon propagator in continuum theory it is impossible to "see" linearly rising potential 
between heavy quarks by lattice QCD simulations [39] not involving some extra (besides 
gluons and quarks) degrees of freedom. This should be considered as a strong lattice evidence 
(though not direct) of the existence and importance of g~^-type excitations of gluon field 
configurations in the QCD vacuum. There also exists direct lattice evidence that the zero 
modes are enhanced in the full gluon propagator indeed [40]. 

e) . Within the distribution theory [30] the structure of the nonperturbative IR singu- 
larities in four-dimensional Euclidean QCD is the same as in two-dimensional QCD, which 
confines quarks at least in the large Nc limit [41]. In this connection, let us note that 

IR singularity is the simplest nonperturbative power singularity in four-dimensional 
QCD as well as IR singularity is the simplest nonperturbative power singularity in 
two-dimensional QCD. The QCD vacuum is much more complicated medium than its two- 
dimensional model, nevertheless, the above-mentioned analogy is promising even in the case 
of the nonperturbative dynamics of light quarks. 

f) . Some classical models of the QCD vacuum also invoke behavior of the gluon fields 
in the IR domain. For example, it appears in the QCD vacuum as a condensation of the 
color-magnetic monopoles (QCD vacuum is a chromomagnetic superconductor) proposed by 
Nambu, Mandelstam and 't Hooft and developed by Nair and Rosenzweig (see Ref. [42] and 
references therein. For recent developments in this model see Di Giacomo's contribution in 
Ref. [14]) as well as in the classical mechanism of the confining medium [43] and in effective 
theory for the QCD vacuum proposed in Ref. [44]. 

g) . It is also required to derive the heavy quark potential within the recently proposed 
exact renormalization group flow equations approach [45]. 

h) . It has been shown in our papers that the singular behavior (2.2) is directly related 
to light quarks confinement and DCSB [27,28]. 

i) . These excitations are also topologically nontrivial since they lead to the nontrivial 
(truly nonperturbative) vacuum energy density which is finite, automatically negative and 
it has no imaginary part (stable vacuum) [21,27,28]. 
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j). Moreover, they nicely saturate the phenomenological values of the topological sus- 
ceptibility, the mass of the t]' meson in the chiral limit and the gluon condensate [25,26]. 

Thus we consider our main Ansatz (2.1) as physically well- motivated. The first term in 
the decomposition (2.1) is the truly nonperturbative part since it vanishes in the perturbative 
limit (/i^ — > 0), when the perturbative phase survives only and simultaneously it correctly 
reproduces the deep IR asymptotics of the full gluon propagator. 



III. THE EFFECTIVE POTENTIAL FOR COMPOSITE OPERATORS 

In this section we are going to analytically investigate the quantum structure of the YM 
vacuum which is due to g~'^-type nonperturbative, topologically nontrivial excitations of the 
gluon field configurations there complemented by the renormalization group improvements 
for the UV region, (2.1). Let us start from the gluon part of the vacuum energy density 
which to leading order (log-loop level ~ h)f\ is given by the effective potential for composite 
operators [19] as follows 



(27r)" 



(3.1) 



where D{q) is the full gluon propagator (1.1) and Do{q) is its free perturbative (tree level) 
counterpart. Here and below the traces over space-time and color group indices are under- 
stood. The effective potential is normalized as V{Dq) = 0, i.e., free perturbative vacuum is 
normalized to zero. In order to evaluate the effective potential (3.1) we use the well-known 
expression, 



Tr\n{DQ^D) = 8 x Indet^D^^D) = 8 x 41n 
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(3.2) 



It becomes zero (in accordance with the above mentioned normalization condition) when 
the full gluon form factor is replaced by its free counterpart. The composition (3.2) does 
not explicitly depend on a gauge choice. 

Going over to four {n = 4) dimensional Euclidean space in (3.1-3.2), on account of (2.1), 
by introducing the appropriate dimensionless momentum variable a.s z = q^ / /i^ and after 
doing some algebra, we finally obtain the following expression for the vacuum energy density. 



:q^ZQ^Ig{zo;b), 



(3.3) 



where 



/ 3 3c 
dz z \ \n z + a + — — ln[3 + z + 
^ 4z Amoz 



3cz 
In bz 



:zo;b) + Fizo;b), (3.4) 



^Next-to-leading and higher terms (two and more vacuum loops) are suppressed by one order of 
magnitude in powers of h at least and are left for consideration elsewhere. 
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and 



Igi^o; b) = -Zq \nbzo - -zl + -Zo - I{zq- h) 

— 3c X X 

Il{z^;h) = -^zl + -^li{bzof + -zlXnXnhzQ + i^Ei{l, -21n6^o), (3.5) 

with a = a + \n.h and a = (3/4) — 2 In 2 while 6 = Infc^o + 3c. In these integrals zq = Qo/fJ''^ 
while b = fi'^/AyM li{x) is the logarithm- integral function. Let us note in advance, 
that here and everywhere below the combination bzo = Qq/^m = is fixed (constant) 
when qQ itself is fixed, i.e., before going to infinity at final stage. Thus, the parameter 
Xq is the truly dimensionless UV cutoff since when the UV cutoff in momentum space 
goes to infinity (^q oo) then it certainly goes to infinity as well. At the same time, 
the variable zq only looks like as the dimensionless UV cutoff. In the limit go ^ cxd it 
might be finite since the nonperturbative scale, /x^, can be affected by the UV divergences 
as well ("bare" parameter). Precisely this takes place within our approach in this case (see 
below). In fact, variable zq determines the perturbative limit {zq oo) at fixed ^q, when the 
perturbative phase survives only in our model (yU 0). Ei{l^ — 21n62;o) is the exponential- 
integral function with unphysical singularity at hzo = 1 as well as the preceding term has. 
In the compositions shown in (3.5), the corresponding integral is defined as follows 

I{zo-h) = j dzzhi[{?, + z)hihz + 2,cz\. (3.6) 



The effective potential at the log-loop level for the ghost degrees of freedom is 

where G{p) is the full ghost propagator and G'o(p) is its free perturbative counterpart. The 

effective potential V{G) is normalized as V{Go) = 0. Evaluating the ghost term egh — 
V{G) in (3.7), we obtain egh = 7r~^kQyQ'^Igh{yo;b), where the integral Igh{yo;b) depends on 
the ghost propagator, which remains arbitrary (unknown) within our approach and t/q = 

In principle, we must sum up all contributions in order to obtain total vacuum energy 
density (the confining quark part of the vacuum energy density is not considered here). 
However, approximating the full gluon propagator by the sum of its deep IR (confinement) 
and UV asymptotics in (2.1), the nonzero terms appear in (3.4) which precisely violates the 
above mentioned normalization condition of the free perturbative vacuum to be zero when 
the parameter b goes to zero {/i — > 0, — oo). In this case the perturbative phase only 
remains within OTir model. These terms arc collected in the expression (3.5) for Ig[zQ]b). 
It also contains two terms which suffer from unphysical singularity - Landau pole - at 
Xq = 1. This should be traced back to (2.1) when — = AyM since we use it in the 
whole range. These terms should be subtracted from the gluon part of the vacuum energy 
density. Using the arbitrariness of the ghost term (3.7), this can be done by imposing the 
following condition, A = Ig{zo;b) + Igh{yo;b) — (recalling that Ay^ = ^'o^o'^ — Qo^o'^)- 
We regularize the gluon contribution to the vacuum energy density by subtracting unwanted 
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terms by means of the ghost contribution, i.e., it plays the role of the corresponding counter 
term. So we define, eg + = e^^^ = ^''^q^z^'^Il^izo-b), on account of A = 0. Thus our 
regularization procedure is in agreement with general physical interpretation of ghosts to 
cancel the effects of unphysical degrees of freedom of the gauge bosons [1,46]. 
Let us introduce the effective potential Clg at fixed go as follows [21,28] 

Qg = Qg{zo; b) = le-^ = ^ x z^'{z', Inbzo - \zl + - 21{z^- 6)}, (3.8) 

where the integral I{zQ\h) is given in (3.6). First of all, it is always useful to factorize 
the scale dependence in the regularized and subtracted vacuum energy density by means 
of the above-introduced effective potential at fixed scale. On the other hand, it makes it 
possible to investigate the structure of the vacuum within our approach since this quantity 
is finite when the UV cutoff goes to infinity (see below) as well as it is free from unphysical 
singularities (see above). These two remarks emphasize the importance of the effective 
potential (3.8). The asymptotic behavior of the effective potential (3.8) depends on the 
asymptotic properties of the integral I{zq; b), (3.6). It is almost obvious that its asymptotics 
at zq 0, oo (and consequently at b 0, oo) to- leading order can be easily evaluated 
analytically by replacement In bz by In bzQ since logarithm is a very slowly varing function. 
Omitting all intermediate analytical calculations, one obtains 

I{zo;b) = IzlHbzo) - \zl + - |^ln{l +~^} + \z',ln{l + |-}, (3.9) 

where c— 1 + (3c/ \nbzo)- Thus in the perturbative limit 2^0 — > oo one has 

n r h\ ^ In bzp - c _^ 

^9izo;b)r.,^^^- — ^^^^^z, , (3.10) 

which clearly shows that in order for the effective potential (3.8) as a function of zq to 
approach zero from below in the perturbative limit, one needs to set 

lnxo = Ac, A>1, (3.11) 

recalling that bzo — Xq by definition. The expression (3.11) is nothing but useful parametriza- 

tion of the UV divergences in terms of the parameter A, which will go to infinity at the final 
stage. At the same time, at zero {zq 0) the effective potential (3.8) behaves as ~ Zq-^. 
Thus as a function of Zq it has a minimum if the condition (3.11) takes place which is the 
case indeed (A goes to infinity at final stage). 

The minimization of the effective potential (3.8) with respect to zq (at fixed bzo), 
dQg{zo; b)/dzo — 0, yields 

z^ + 4l2{zo; b) = ^zo + 2z^o Hi^ + ^0) In 6^0 + Sc^o]- (3.12) 

It is easy to show that to- leading order in the A ^ 00 limit, it can be analytically evaluated 
as follows: 

(A - 1)(A + 3)zo = AXHn{l + ^^^|^}. (3.13) 
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In the derivation of this expression we have aheady used expression for c (see text after Eq. 
(3.9)) and (3.11). Its solution in the A — oo hmit is 

;,-n(A ^ oo) = {^} . (A ^ oo) = 2.2, (3.14) 

^ > mm 

i. e. the point 2;™" = 2.2 is the point of accumulation of minimums in the UV limit of the 
effective potential (3.8). This effect shows that the parameter h goes to infinity in the same 
way as the UV cut-off xq in order for their ratio to go to the above displayed finite limit. 
This means that initial mass scale parameter is a "bar" one. It depends on the UV cut-off 
A, i.e., /i = /i(A) and consequently h = &(A)|. In what follows we will introduce the UV 
renormalized fi and consequently the UV renormalized b. To-leading order in the A ^ 00 
limit the effective potential (3.8), on account of the "stationary" condition (3.13), becomes 

«.(^.^») = ^{j^o-'-Hl + ^^)}. (3.15) 



IV. THE YM VACUUM STRUCTURE 

Evaluating the effective potential fig(A, 2;™"(A)) (3.15) numerically in the vicinity of the 
point of accumulation of the minimums, z™"(A 00) 2.2 (see (3.14)), one obtains 
^^(A —>■ 00) = —0.0263, i.e., it is finite in this limit. Recalling that e^*^^ = e^^^(A), then from 
Eq.(3.8) it follows 

e7(A) = -0.0263go'(A). (4.1) 

On the other hand, 

q',{X) = xl{X)AU, = e'^^A'yM, (4.2) 

by definition and because of (3.11). Thus the vacuum energy becomes 

e^'^9(A) = -0.0263A^A^ x e^^' (4.3) 

and it is badly divergent in the UV limit (A 00) because of the exponential A-factor. This 
result should be expected since we have explicitly kept (did not subtract) the UV peace in 
Eq. (2.1) (for general discussion see Ref. [21]). 

However, despite badly divergent vacuum energy density, our method makes it possible 
to establish a finite relation between two scale parameters of our model. Indeed, as was 
mentioned above, the parameter b diverges exactly as the UV cutoff xq = e^'^ (see Eq. 



^Let us note that it does not make any sense to introduce finite (from the very beginning) the 
nonperturbative mass scale parameter. The problem is that in any case it will be affected by the 
UV divergences since we have explicitly kept the UV term in Eq. (2.1) (did not subtract it). At the 
same time, the asymptotic scale parameter, Ay^f, is finite (from the very beginning) by definition. 
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(3.11)) itself in order to get final result in the A — > oo limit. So let us introduce now the UV 
finite b and /i as 

b = 6(A)e-^" (4.4) 

and 

A^P = //(A)e-^^/^ (4.5) 
respectively. Then from Eq. (3.14) one obtains 



Ac 

Zr, 



= 2.2 (4.6) 



6(A) 66^^= 

in the A — > oo limit. On the other hand, by definition 

in the same hmit. This finally leads to 

Aym = 1.48324 X Knp- (4.8) 

A few remarks are in order. Let us underline that the above- displayed numerical result 
(4.8) was obtained in a gauge invariant way, i.e., this number did not explicitly depend on a 
gauge choice. Also this result will not be changed in the next-to-leading order approxima- 
tion to the running coupling constant in the UV limit since they were obtained to-leading 
dominant order limit. Evidently, it should not, of course, depend on how one introduces the 
UV cutoff by making use of different scale parameters. For example, it can be defined as 
Xq = QqI Kym (from the very beginning) while the parameter b remains the same. Then the 
minimization of the effective potential VLg{xQ]b) = e^^^(a;o; 6)Ay^ with respect to b finally 
yields the same numerical relation between the nonperturbative and asymptotic scales, (4.8) 
as it should be in principle (for details see Ref. [47]). It is worth noting that the minimization 
with respect to xq itself in this case leads to trivial zero (see Ref. [21]). 



V. DISCUSSION 

The main problem now is the numerical value of the YM asymptotic scale parameter, 
Aym = A (for brevity and further aims), since it can not be experimentally determined 
in the YM theory alone. So its numerical value is to be taken from lattice approach or 
phenomenology. Unfortunately, lattice estimates of the QCD flavorless asymtotic scale pa- 
rameter strongly depend on the renormalization scheme chosen for calculation. In Ref. [48] 
in the MS scheme it is expressed in terms of Sommer's scale Tq, namely A;^ = 0.636(54)/ro. 
Than Eq. (4.8) becomes A^p — 0.4288r^^. In physical units, using tq = 0.5 /m, one nu- 
merically gets, Ajvp = 171.5 MeV. However, as it was emphasized in the above mentioned 
Ref. [48] , the conversion to physical units is rather ambiguous in the pure gauge theory. Ap- 
proximately the same numerical result for Aj^ in terms of the string tension was previously 
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obtained in Ref. [49]. In Ref. [50] the YM scale parameter was calculated first in the MOM 
scheme and then it was matched to the MS scheme. The corresponding numerical results 
(here and everywhere below up to central values, for simplicity) arc: ^mom ~ 0-88 GeV and 
^TTs ~ '-'•31 GeV . Substituting these values into the Eq. (4.8), one arrives at the following 
estimates, K^^^ = 593.3 MeV and = 210 MeV . The most recent lattice calcula- 
tion of the flavorless QCD asymptotic scale parameter is being presented in Ref. [51] are: 
A^o^ = 0.361(6) GeV and ^M^M = 0.349(6) GeV while A^ = 0.412 GeV. Substituting 
this result into the Eq. (4.8) for example for the last value, one gets, Ajlfp = 277.7 MeV . 
In Ref. [52] the vacuum energy density of the covariant model of Euclidean QCD with 
= 3 quark flavors on a finite Z^- dimensional symmetric hypcrtorus was calculated. 
Adjusting the recently calculated value of the gluon condensate by Narison [53] and the 
numerical value of ag from r decay, in the above-mentioned paper [50] it was estimated 
that A^ ~ 500 — 600 MeV (compatible with the standard value of the gluon condensate 
[54]). From the renormalization group equations, however, it is well-known that the nu- 
merical value of the flavorless QCD asymptotic scale parameter should be bigger than its 
full QCD counterparts, for example, the above estimated A^^. Thus roughly estimating 
Aym ~ 600 - 700 MeV in Eq. (4.8), one obtains Anp ~ 404.5 - 472 MeV. Contrary to the 
above-discussed cases, in Refs. [27,28] the nonperturbative scale parameter was calculated 
itself in chiral QCD within our approach. Because of the scale-setting scheme in order to 
relate it to a good physical observeble from full QCD , its numerical value depends on a 
bounds for the pion decay constant in the chiral hmit. The corresponding bounds for the 
nonperturbative scale parameter, we have obtained, were: 540.6 < fl < 578.4 MeV. Then 
from Eq. (4.8) it immediately follows 801.8 < Aym < 858 MeV by identifying p, with 
Anp for simplicity. Here we are observing the standard situation when estimates from phe- 
nomenology (continuum theory) are usually bigger than lattice estimates of the same scale 
quantities. 

In any case, our solid result is the relation (4.8) which is to be better understood as 
the ratio between two scale parameters from lattice simulation point of view. In this way, 
in principle, it should be free from lattice artifacts since it does not require the conversion 
to physical unuts which still remains rather ambiguous in the pure gauge theory as it was 
mentioned above. 

VI. CONCLUSIONS 

In summary, we have formulated new quantum model of the QCD vacuum using the 
effective potential approach for composite operators. It is based on the existence and im- 
portance of such kind of the nonperturbative, topologically nontrivial excitations of gluon 
field configurations (due to self-interaction of massless gluons only, i.e., without explicit in- 
volving some extra degrees of freedom) which can be effectively correctly described by the 
g~^-type behaviour of the full gluon propagator in the deep IR domain. The UV part of the 
full gluon propagator was approximated by asymptotic freedom to-leading order log term of 
the running coupling constant. By explicit retaining the UV tail in the decomposition of the 
full gluon propagator (2.1), the vacuum energy density remains badly divergent, of course 
(for general discussion see again Refs. [21,22] and Shuryaks's review in Ref. [55]). However, 
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despite this, we have formulated a method to estabhsh a finite (in the UV hmit) relation 
(4.8) between two scale parameters of our model. We propose to minimize the properly reg- 
ularized and subtracted effective potential at a fixed scale (3.8) as a function of the variable 
Zq. It has a clear physical meaning when it goes to infinity then the perturbative phase 
survives only within our approach to nonperturbative QCD. This is possible to do because 
the above-mentioned effective potential is finite in the UV limit within our model, i.e., this 
is a specific property of our Ansatz for the full gluon propagator (2.1). Equivalently, one 
can minimize the auxiliar (again properly regularized and subtracted) effective potential as 
a function of the parameter which has a clear physical meaning. When it is zero then again 
the perturbative phase survives only. 

For the first time, we have expressed the asymptotic scale parameter as pure number 
times the nonperturbative scale, namely (4.8). We don't know any other YM quan- 
tum/classical vacuum model which possesses the same property. This can be considered 
as one more argument in favor of our model. In principle, in QCD one of the final goals is 
precisely to express all scale parameters (for example, hadron masses) as pure numbers times 
some characteristic scale of the theory. Despite the perturbative vacuum cannot be the QCD 
true ground state [56], nevertheless the existence of such kind of relation is a manifestation 
that "the problems encountered in perturbation theory are not mere mathematical artifacts 
but rather signify deep properties of the full theory" [57]. In support of this statement, let 
us underline one more that the running coupling with scaling violations not only in the UV 
but in the IR as well is the most important feature of QCD which dratically differs this 
theory from QED. The former is certainly related to asymptotic freedom while the latter 
should be apparently related to confinement, in particular quark confinement. Behind both 
types of scaling violation (in general, with their own scale parameters) is the self-interaction 
of massless gluons only within our approach to nonperturbative QCD in general and to its 
ground state in particular. Thus it suggests a unified description of both phenomena, so 
such kind of relation between two scale parameters, in principle, should be expected. At the 
same time, it does not mean that they are dependent from each other. It means rather that 
there exists the procedure (method) how to establish connection between them. Precisely 
this point we wanted to emphasize in the Introduction. 
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